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We describe a strategy for a non-perturbative computation of the b-quark mass to leading order in 1/m in 
the Heavy Quark Effective Theory (HQET). The approach avoids the perturbative subtraction of power law 
divergencies, and the continuum limit may be taken. First numerical results in the quenched approximation 
demonstrate the potential of the method with a preliminary result m^^(4GeV) = 4.56(2)(7) GeV. In principle, 
the idea may also be applied to the matching of composite operators or the computation of 1/m corrections in 
HQET. 



1. THE PROBLEM 



2. STRATEGY 



Since the b-quark is heavier than presently 
achievable inverse lattice spacings, 1/a, effective 
theories (HQET, NRQCD, ...) are used for its 
numerical treatment. These theories are based 
on some large mass expansion and hence do not 
have a chiral symmetry to protect the quark mass 
from additive renormalization. The relation 



Z(m{^"<= + 6m) = mb 



(1) 



between the bare quark mass mj^"'' and a renor- 
malized b-quark mass, mb, therefore contains an 
additive term, Sm. For dimensional reasons it is 
linearly divergent. 
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A perturbative approximation to f{go) results in 
a truncation error in eq. (Q) which diverges in 
the (7o ~* limit. The continuum limit can- 
not be taken. Available results for the b-quark 
mass in the static approximation [Q are limited 
by this fact and consequently are obtained at lat- 
tice spacing a w 0.1 fm. 

While the numerics may be improved by 
the computation of higher order fi a non- 
perturbative strategy to compute or avoid dm is 
needed to solve the problem. 



We treat the b-quark field in static approxima- 
tion with action 5*11 = i'h^oi'h (note that Sm 
is not included in the action and consult for 
any unexplained notation) and consider a times- 
lice correlation function Cstat(a;o) of some inter- 
polating field for the B-meson. Its large time de- 
cay Cstat(a;o) ~ i?exp(— i?xo) defines the static 
bare binding energy E, related to the B-meson 
mass, rriB: via 



mB^E + ml''''' -}-0(l/TOb). 



(4) 



♦Presented by R.S. at "Lattice 2001" 



2.1. Matching 

In order to replace mb^''° by the renormalized 
quark mass we consider another such relation, 
which we may take to be the condition match- 
ing the QCD quark mass to the one of HQET. It 
reads 

r,ei(L,M,5o) = r,tat(i,ffo) + m|^^''=(M,go) (5) 

with +0(1/^^) 
rstat(i,5o) = -^(5o + 9j)ln[Cstat(a;o)] _l (6) 

xq— 2 

and Trc\{L, Af, go) defined in the same way but for 
a relativistic quark with renormalization group 
invariant (RGI) quark mass, M . Its relation to 
the bare PCAC mass and the bare mass in the 
Lagrangian of the 0(a)-improved theory is known 
non-perturbatively At this stage L may be 

any finite length scale. 



In principle the function m{^"°(M, go) can be 
obtained by evaluating eq. (^) for some range of 
M. It may then be inserted into eq. and solved 
for M to obtain the RGI b-quark mass Mb. 

In practice, not much has been achieved yet 
because an implementation of eq. (^) is not pos- 
sible: the small lattice spacings necessary for the 
relativistic theory are not available. 

2.2. The use of finite volume 

The important idea is to consider C(xo) to be 
a correlation function in finite volume of linear 
extent L, which may assume several values L = 
Li. We then have (suppressing the argument go) 

me - £;-r,tat(io)+rrei(Lo,Mb) (7) 
= Ai;a + A£;b + r,ei(Lo,Mb), (8) 

ASa = £;-rstat(i«), (9) 

A£;b = r3tat(i„)-rstat(io). (lo) 

In both energy differences AiJ^, defined in the 
static theory, the unknown m^^^'^ (and thereby 
(5™, too) cancels, and the continuum limit of 



r2a — —LqAEs, , fib — —LqAEi, 



(11) 



exists. Choosing furthermore Lq such that 
Lo-^b 3> 1 and MbO ^ 1 can be achieved at 
the same time, the function 



^01(2) = iorrci(io, M) , Z = LoM , 



(12) 



may be evaluated (in the relativistic theory) and 
extrapolated to the continuum limit. Then Afb 
can be determined by solving 

aci(ioMb) = f^a + f^b + f^c , (13) 

with the experimental input (spin averaged mass) 

fic = LoniB, =LoX 5405 MeV. (14) 

On the one hand, to treat the relativistic b-quark, 
Lq may not be too large and on the other hand, 
to be able to compute Ai?a, Ln may not be too 
small. To bridge this gap, a step scaling function 



ar(u) = lim Er(u,a/L), 



(15) 
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Figure 1. Step scaling functions for u 
2.4484, 2.77, 3.48 and (for ar only) u = 2.1 . 



depending on a coupling, g'^{L), renormalized at 
scale L, is introduced in the spirit of [Q. For 
Li = 2*Lo, Ui = g'^{Li), one then has 



5]2-^-Vr(uO. 



(16) 



j=0 



We now demonstrate the accessibility of the con- 
tinuum limit with preliminary 

3. RESULTS 

in the quenched approximation. We choose 
Schrodinger functional (SF) b.c.'s and specifically 



C{xo) = fAixo) 



(17) 



Sr(^/,f) = 2L[rstat(2L)-r3tat(i)]„=g2(i) 



9 — 0.5, T ~ L and vanishing background field; 
see for the definition of /a in the static ap- 
proximation and with a relativistic b-quark. 
is taken to be the standard SF coupling [^J|]; 
we set Lq = 1.436ro/4 « 0.2 fm and the light 
quark mass to zero, except in the large volume B- 
mcson correlation function where, in accordance 
with eq. (p^, it is set to the strange quark mass 
0. All quantities are 0(a)-improved. 

After continuum extrapolation of Sr and a 
slight interpolation of the function a{u) (see 
Fig. ^ we estimate {n = 2) 

Ob = 0.10(1). (18) 

The binding energy E is computed in a (1.5 fm)^ x 
2.3 fm box with SF b.c.'s and a new technique to 
suppress excited states 1^. From the plateau in 
Fig. H we obtain 

r^a = -0.369(6) at /? = 6/5^ = 6.0 . (19) 
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Figure 2. Plateau for fJa with E = Ecs{xo) = 
a-i[ln{C(a;o)) - ln(C(xo - a))] and a « 0.1 fm. 
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Figure 3. riroi and graphical solution of eq. 



Finally, the missing piece of the puzzle, r2i.oi(z), 
is plotted for various L/a- and z- values in Fig. |[ 
where also the solution of eq. ( p^ ) is illustrated. 
Setting To = 0.5 fm, we end up with 

Mb = 7.01(3)(10)GeV, (20) 
m^(4GeV) = 4.56(2)(7) GeV , (21) 

where the second, presently dominating uncer- 
tainty originates from the renormalization factors 
and improvement coefficients needed to fix M in 
the relativistic theory. 

4. COMMENTS, GENERALIZATIONS 

Our computation of is valid up to (rel- 
ative) 1/Mb corrections. Generically these are 
of the form fi/M^, with /i « 0.5 GeV, a typical 
QCD scale. However, also explicitly introduced 
external scales may enter. In our strategy only 
the scale Lq where the matching is performed 
appears. We therefore do not choose Lq much 
smaller but keep 1 /Lq « 1 GeV. 

Since fia, eq. (pj|), has only been evaluated at 
one value of a, our present estimate does not yet 
reflect a final result in the full continuum limit. 
Nevertheless, it is already useful to compare this 
with earlier results which were obtained at similar 
lattice spacings but perturbative Sm. Within the 
present errors we agree with them |l] . Our results 
still need a more detailed analysis of discretiza- 
tion effects also in f2b,^2rci- The uncertainties 
of improvement coefficients and renormalization 
constants in the range of 170 relevant for Fig. ^ 
are presently contributing the dominant 70 MeV 
error in eq. (pl|). It can and should be reduced. 

A relevant issue not mentioned in the litera- 
ture of lattice determinations of Wb is the effect 



of the scale ambiguity in the quenched approxi- 
mation (see e.g. sect. 6 in 0). We found that a 
10% ambiguity in the scale roughly corresponds 
to a 100 MeV error in the quark mass. 

A theoretically interesting point is that in our 
strategy power law divergences in the effective 
theory are removed by a non-perturbative match- 
ing to relativistic QCD in finite volume. The 
same strategy may e.g. be applied to the renor- 
malization of the static axial current and to 1 /Mb 
corrections in HQET. Numerical experiments are 
needed to test the practicability of such general- 
izations, which open up exciting possibilities for 
applications of HQET without perturbative (in 
the QCD-coupling) truncation errors. 
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